o 

vO 

CO 

I 

iT) 

vO 

I 

cx: 

H 

I 

UJ 


I 

>•  I 
a.  I 

8  I 
a  I 

i-f  I 
tL,  I 

I 

H  I 
H  ( 
(fi  I 
uj  I 


ESD-TR-65-369 


ESD  ACCESSION  LIST 

ESTI  Call  No._  AL  47351 

Copy  No.  /  0,  t  cys. 


TM-03967 


THE  REDUCTION  OF  THE  VARIANCE  BY 
LEAST  SQUARES  POLYNOMIAL  APPROXIMATION 


i 

i 


TECHNICAL  REPORT  NO.  ESD-TR-65-369 
SEPTEMBER  1965 


H.  C.  Joksch 


ESD  RECORD  COPY 

RETURN  TO 

SCIENTIFIC  &  TECHNICAL  INFORMATION  DIVISION 
(ESTI),  BUILDING  1211 


Prepared  for 

DIRECTORATE  OF  COMPUTERS 
ELECTRONIC  SYSTEMS  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 

L.  G.  Hanscom  Field,  Bedford,  Massachusetts 


Project  7070 
Prepared  by 


THE  MITRE  CORPORATION 
Bedford,  Massachusetts 
Contract  AF  19(628)-2390 


Copies  available  at  Clearing  House  for  Federal 
Scientific  and  Technical  Information  (formerly  Office 
of  Technical  Services). 


Qualified  requesters  may  obtain  copies  from  DDC. 
Orders  will  be  expedited  if  placed  through  the  librarian 
or  other  person  designated  to  request  documents 
from  DDC. 


When  US  Government  drawings,  specifications,  or 
other  data  are  used  for  any  purpose  other  than  a 
definitely  related  government  procurement  oper¬ 
ation,  the  government  thereby  incurs  no  responsi¬ 
bility  nor  any  obligation  whatsoever;  and  the  fact 
that  the  government  may  have  formulated,  fur¬ 
nished,  or  in  any  way  supplied  the  said  drawings, 
specifications,  or  other  data  is  not  to  be  regarded 
by  implication  or  otherwise,  as  in  any  manner 
licensing  the  holder  or  any  other  person  or  corpo¬ 
ration,  or  conveying  any  rights  or  permission  to 
manufacture,  use,  or  sell  any  patented  invention 
that  may  in  any  way  be  related  thereto. 


Do  not 


return 


this  copy. 


Retain  or  destroy. 


ESD-TR-65-369 


TM-03967 


THE  REDUCTION  OF  THE  VARIANCE  BY 
LEAST  SQUARES  POLYNOMIAL  APPROXIMATION 


TECHNICAL  REPORT  NO.  ESD-TR-65-369 
SEPTEMBER  1965 


H.  C.  Joksch 


Prepared  for 

DIRECTORATE  OF  COMPUTERS 
ELECTRONIC  SYSTEMS  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 

L.  G.  Hanscom  Field,  Bedford,  Massachusetts 


Project  7070 
Prepared  by 

THE  MITRE  CORPORATION 
Bedford,  Massachusetts 
Contract  AF  19(628)-2390 


b 


ABSTRACT 


Least  squares  approximations  obtain  parameters  with  a  variance 
lower  than  those  of  the  data  from  which  they  are  obtained.  A 
least  squares  polynomial  approximation  to  observations  may  be 
used  to  obtain  '^smoothed"  values  of  the  observations  or  to  make 
predictions.  The  reduction  of  the  variance  achieved  by  this 
process  is  determined  for  several  special  cases.  Some  properties 
of  the  Gram -polynomials  necessary  for  the  analysis  are  derived. 
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SECTION  I 
INTRODUCTION 


APPROXIMATION  TECHNIQUES 

If  a  function  is  empirically  given,  e.  g.  ,  by  discrete  observations,  it  is  frequently 

0  desirable  to  represent  these  observations  by  analytical  expression,  e.  g.  ,  polynomials.  If 

the  values  of  the  function  are  not  exactly  known,  but  subject  to  error,  then  it  is  not  reason¬ 
able  to  obtain  an  exact  representation.  It  is  more  appropriate  to  use  an  approximate  repre¬ 
sentation,  where  the  differences  between  the  observations  and  the  approximation  are  so 
small  that  they  can  be  considered  as  errors  of  the  observations  rather  than  as  inaccuracies 
of  the  representation.  If  the  ’true"  nature  of  the  empirical  function,  e.  g.  ,  the  degree  of  a 
polynomial,  is  known,  it  can  be  used  to  approximate  the  observation  and  determine  its 
parameters.  The  differences  between  approximation  and  observation  can  be  considered  as 
pure  errors  of  the  observations. 

One  way  to  find  the  "best"  approximation,  which  minimizes  an  overall  measure  for  the 
error  of  the  approximation,  is  a  least  squares  approximation.  The  randomness  of  the 
errors,  however,  precludes  an  approximation  method  that  will  find  the  "true"  function  and 
separate  it  from  the  observational  errors.  Thus,  some  part  of  the  observational  errors 

*  can  always  influence  the  approximating  function,  and,  consequently,  the  differences  will  con¬ 
tain  not  only  observational  errors  but  also  inaccuracies  of  the  approximation.  The  larger 

•  the  number  of  errors,  the  more  likely  the  probability  that  the  random  effects  will  cancel  out 
and  that  the  approximating  function  will  be  close  to  the  "true"  function.  On  the  other  hand, 
the  more  the  parameters  necessary  to  describe  the  approximating  function,  the  less  the 
reduction  of  this  random  error  will  be.  These  functional  dependencies  are  the  subject  of 
this  analysis. 

LEAST  SQUARES  APPROXIMATIONS 

This  analysis  is  restricted  to  an  important  special  case:  the  values  of  the  empirical 
function  ("observations")  y^  are  given  for  n  equidistant  values  of  the  argument 
X,  X  =  0  , . .  n-1,  and  the  approximating  function  is  a  polynomial  y(x). 


$ 
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Some  important  applications  of  least  squares  fitting  include: 

(a)  smoothing  a  series  of  values; 

(b)  estimating  more  precisely  the  present  position  of  an  object  during  tracking;  and 

(c)  determining  the  expected  position  y(x)  of  the  object  being  tracked  and  its  vari¬ 
ance  at  the  time  of  the  next  observation. 

In  (a),  the  polynomial  values  at  the  midpoint  of  the  components  are  used  as  ’'smoothed^' 
values;  in  (b) ,  a  polynomial  of  appropriate  degree  is  fitted  to  the  last  observations,  and 
y(x)  at  the  endpoint  is  used  as  the  last  position.  Thus,  exact  expressions  and  bounds  for 
the  variance  of  the  polynomial  can  be  given  for  (a) ,  (b) ,  and  (c)  at  the  midpoint,  the  end¬ 
point,  and  at  the  first  equidistant  point  outside  the  given  interval,  respectively. 

General  expressions  for  the  variance  are  well  known,  but  their  use  requires 

extensive  numerical  calculations.  The  dependency  of  the  variance  on  the  number  of  observa¬ 
tions  and  on  the  degree  of  the  approximating  polynomial  is  not  obvious.  Special  studies, 

[1]  [3  4  5] 

sometimes  incidental  to  other  problems,  have  been  made  by  Cowden,  Guest,  *  ’ 

Proschan,  and  Smith.  Smithes  problem  is  slightly  different;  she  assumes  that  the 

observations  are  uniformly  spread  over  the  whole  interval.  Her  exact  results  for  the  mid- 

and  endpoints  are,  therefore,  only  approximations  for  the  problem  discussed  in  this  report. 
[4] 

Guest  obtains  the  same  approximations  in  a  much  easier  way  by  the  use  of  Legendre 
polynomials.  Proschan  finds  the  same  approximation  for  the  endpoint  by  use  of  determi¬ 
nants,  similar  to  Smith.  However,  these  approximations  are  good  for  large  n  only;  for 

[3  5 1 

small  n,  they  can  be  very  poor.  Guest  ’  and  Cowden  give  tables  for  the  determina¬ 
tion  of  the  variance  of  y(x) ,  based  on  numerical  calculations.  Using  Guest ^s  tables, 
one  has  to  notice  that  k  =  1  is  not  the  endpoint  of  the  interval,  but,  rather,  k  =  (n-l)/n. 
Some  checks  showed  good  agreement  with  the  results  obtained  from  our  exact  formulas. 


Numbers  in  brackets  refer  to  citations  in  the  Bibliography. 
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SECTION  n 

APPROXIMATION  BY  ORTHOGONAL  POLYNOMIALS 


[5, 6, 7,8] 


Let  n  values  y 
X  =  0  , . .  n-1.  Let  y(x) 
that 


(’’observations’')  be  given  for  equidistant  arguments  (’’times”) 
be  a  polynomial  of  degree  m  to  be  fitted  to  the  y^  such 


4 


n-l 
x=0 


2 

=  mm. 


(1) 


If  the  y^  are  subject  to  random  variations,  then  the  coefficients  of  the  polynomial  are  also 
subject  to  those  variations  and,  consequently,  the  value  of  y(x)  for  any  given  x. 


Least  squares  approximations  are  simplified  by  the  use  of  orthogonal  polynomials 
(l>^(x).  They  are  characterized  by  the  orthogonality  relations 


I 


W(X)  (X)  (X)  = 


if  r  =  s 
if  r  ^  s 


(2) 


In  the  case  of  equal  weights,  w(x)  =  1,  and  equidistant  arguments,  x  =  0  . . .  n-l,  they 
determine  the  Gram -polynomials  (sometimes  called  Chebysheff -polynomials;  but  this  name 
is  more  commonly  used  for  other  orthogonal  polynomials).  An  explicit  expression  for 
them  is 


i=0 

where  x  =  x(x-l)  . . .  (x-i+1) ,  and  the  sum  of  their  squares  is 

n  n+1  n+2  n+r 

\  2r+l  *  n-l  *  n-2  *  ’  *  n-r 
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The  approximation  which  satisfies  Equation  (  1  )  is  given  by 


where 


(  5 


Y 

r 


From  (  5  )  and  (  6  )  we  obtain 


y^c(x,5)  , 


(6 


(  7 


if  we  define 


C(x,§)  = 


(8 


Equation  (  7  )  is  very  convenient  for  practical  use,  since  it  is  linear  in  the  y 


numerical  values  of  the  c(x,§),  for  several  n  and  m,  are  given  by  Cowden, 
Hildebrand,  and  Milne.^^^ 


The 

[1] 
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SECTION  m 

A  GENERAL  PROPERTY  OF  THE  REDUCTION  FACTOR 


» 


« 


If  the 

the  variance  of 


are  uncorrelated  random  variables,  all  with  the  same  variance  a, 

[5  7 1 

y(x)  is  known  ^  to  be 


then 


V(x)  =  a 


(9) 


The  term  R(x)  -  V(x)  /a  may  be  called  the  reduction  factor.  It  can  be  calculated  for  any 
X  from  Equation  (  9  )  as,  e.g.  ,  in  Cowden  and  Guest.  However,  if  one  uses 
Equation  (  7  )  to  obtain  the  polynomial  value  y(x)  for  any  of  the  given  equidistant  arguments 
X,  it  can  be  obtained  in  a  simpler  way.  Comparison  of  Equations  (  9  )  and  (  8  )  shows 


R(x)  =  c(x,x)  .  (  10  ) 

Verbally  expressed:  the  reduction  factor  for  the  variance,  achieved  by  a  least  squares 
polynomial  approximation  to  equally  weighted  equidistant  observations,  at  any  of  the  observa¬ 

tion  points,  equals  the  coefficient  of  the  observation  at  this  point  in  the  linear  combination  of 
the  observed  values  giving  the  polynomial  value  at  the  point  under  consideration.  This 

gives  immediately,  without  calculations,  the  reduction  factor  if  one  uses  the  numerical 
expression  for  (  7  )  as  given  in  Cowden,  Hildebrand,  and  Milne. 

For  a  special  case,  the  midpoint  of  an  uneven  number  of  observations,  this  result  has 

rsi 

been  obtained  by  Milne. 


% 
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SECTION  IV 

THE  REDUCTION  FACTOR  AT  THE  MIDPOINT 

The  midpoint  of  the  interval  x  =  0  . ,  .  n-1  is  x  =  (n-1)  /  2.  The  value  of  <j)^  at 
this  point  is  derived  in  the  Appendix.  We  substitute  Equations  (  4  )  and  (  46  )  into  (  9  ) , 
write  r  =  2s,  since,  in  effect,  we  have  to  sum  the  even  terms  only,  and  obtain 


— 1  -  1- 

3 

2s-l| 

2 

fn+1 

n+3 

n+2s-l' 

^  4s+l 

n-1 

n-2s 

2  2  • 
/  s=0  ' 

4  •  • 

2s  1 

[n-2  * 

n-4  ■  ■ 

n-2s 

n 

n+1  ‘ 

n+2s 

where  t  is  the  largest  integer  such  that  2t  <  m.  Since  R  =  R  ,  we  restrict  our 
further  arguments  to  m  =  2t.  Extensively  written.  Equation  (  11  )  becomes 


fn-ll 

|_  1 

1+5 

^l| 

r^ 

2t 

2  J 

1  n 

n 

i2  2  . 

n  -1 
2 

n  -n 


+  9 


1  3 

2  •  4 


'  0^-4 


n^-16 


H  ("• 

3 

\2  2 

2t-l  n  -1 

1 

CD 

2  2 
n  -(2t-l) 

4  *  * 

’  2t  2  ’ 

2  “ 

*  2  .^2 

'  n  -4 

n  -16 

n  -  (2t) 

(12) 


This  is  an  exact  expression  for  the  reduction  factor.  To  obtain  a  simple  estimate  for  it, 
we  define 


S 


2t 


+ 

4t+l 

>  1 

1  1 

2t-l’ 

2  ’  4  ” 

2t 

(13) 


We  have,  obviously, 


—  <  R 
n  2t 


n-1 


®2t  n^-1  n^-9 


2  2 
n  -4  n  -16 


2  2 

n  -(2t-l) 

2  2 
n  -(2t) 


(14) 
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It  can  be  shown  by  induction  that 


s  -  ll  1  2t+l]  ^ 

°2t  ^  [2  ■  4  ■  ■  ■  2t  I 


A  comparison  with  Wallis’  product, 


1  =  Lim  1  3  3  5  2^  2t+l 

TT  t=oo  2  '  2  ‘  4  *  4  “ '  2t  ‘  2t 


(15) 


gives  the  estimate 


^(2t+l)  <  <  I  (2t+l) 


(16) 


Combining  (  14  )  and  (  16  ) ,  we  have 


2 

TT 


2t+l 


<  R 


2t 


n-1 


2  2  2  2 
3  2t+l  n  -1  n  -9  n  -(2t-l) 

4*n’2  *2  *'*2  2 

n  -4  n  -16  n  -(2t) 


(17) 


For  large  n,  the  terms  containing  n  are  very  nearly  1,  and,  since  2/tt  =  0.  64,  the 
range  for  given  by  Equation  (  17  )  is  fairly  small.  However,  if  it  is  not  small  com¬ 

pared  to  n,  an  approximation  of  Equation  (  12  )  may  be  worthwhile.  We  notice  that  the 
factors 


4t+l 


1  1  2t-l 

2  ’  4  ‘  “  2t 


in  Equation  (  12  )  approach  4/n  =  1.  27  very  rapidly.  The  first  of  them,  5(1/2)  =  1. 25, 
is  already  very  close  to  the  limit.  Therefore, 


nR, 


2t 


n-l 

”2  2 

n  -1  ,  n  -1 

n^-9  ^ 

n^-1 

n^-(2t-l)  ^ 

2 

2  ^  2  ■ 
n  -4  n  -4 

11^-16 

11^-4 

-  (2t)  ^ 

(18) 
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If  we  expand  the  terms  in  the  bracket  and  neglect  all  terms  of  higher  order  than  n, 


we  obtain 


nR 


2t 


t(t+l)  (4t+5) 
6n^ 


(19) 


Guest  ^"^^and  Smith 
m  is  comparable  to  n. 
and  approximate 


obtain  the  estimate  R^^  ^^2t 


n.  It  is  quite  good,  except  when 
For  a  first  estimate,  one  may  even  make  use  of  Equation  (  16  ) 


.7  (2t+l)  jn. 


COROLLARY 

th 

If  a  polynomial  of  (n-1)  degree  is  fitted  to  n  observations,  an  exact  fit  is 
achieved,  and,  therefore,  R  =  1  for  the  arguments  of  the  observations.  For  odd  n, 
the  midpoint  (n-1)  /2  is  one  of  the  given  points.  Therefore,  we  have  the  identity 


# 


n  =  H-5  - 


2  n  -1 
0^-4 


+  9 


1  3 

2  ‘  4 


2  2 
'2  n  -1  n  -9 


n^-4  n^-16 


+  .  , 


+  (2n-l) 


1  3 

2  ■  4 


•  2  2 

n-2  2  n  -1  n  -9 


n-1 


n^-4  n^-16 


4n-4 

2n-l 


(20) 


for  all  odd  n. 
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SECTION  V 

THE  REDUCTION  FACTOR  AT  THE  ENDPOINT 


At  the  endpoint  of  the  interval,  x  =  0  (the  other,  x  =  n-1 ,  gives  the  same  results) , 
we  have  ^  ^  =  +1 ;  therefore ,  Equation  (  9  )  combined  with  Equation  (  4  )  immediately 
gives 


=n 

m  n 


,  „  n-1  ^  n-1  n-2  n-1  n-m 

1+3  — -  +  5  — T  •  — +  . . .  +  (2m+l)  ’  •  . .  •  • 

n+1  n+1  n+2  '  n+1  n+m 


(21) 


If  we  consider  that  1+3+5  . . .  +(2m+l)  =  (m+1)  ,  we  have,  obviously. 


n  m  ^  ^  n  *  n+1  *  n+2  *  *  *  n+m 


(22) 


is  always  one  of  the  given  points.  Therefore,  R^  (0)  <1  is  obvious  from  Equation  (  21  ). 


The  argument  used  in  the  corollary  of  SECTION  IV  gives  R^  (0)  =1,  since  the  endpoint 
is  always  one  of  the  given  points.  Thereto 
This,  combined  with  Equation  (  22  ) ,  gives 

(m+1)  “  n-1  n-2  n-m 


(m+l)  I  _  (m+l)  n-1  n-2 

Min  ( 1 ,  )  >  R  (0)  >  i ^  .  - -  .  — "  . . . 

\  n  ( —  m '  '  n  m+1  n+2 


n+m 


(23) 


For  small  n,  these  bounds  can  be  quite  far  apart.  To  obtain  an  approximation,  we  expand 
the  fractions  in  the  bracket  in  Equation  (  21  )  in  series  and  omit  terms  containing  n  and 
higher  powers.  This  gives 


(m+1)^  -  ~  j^3  •  1  +  5  (1+2)  +  . . .  +  (2m+^(l+2+  . . .  +  m)  j  . 


(24) 


The  bracket  can  easily  be  evaluated,  and  we  obtain,  by  combining  Equations  (  21  )  and  (  24  ), 

.2  r 


m  n 


I  (m+2) 
2n 


(25) 
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[4]  [9]  [lO]  2/2 

Guest,  Prose han,  and  Smith  ^  give  the  approximation,  (0)  ^^^(m+l)  j x\. 

This  approximation  is  good  for  very  large  n  only.  It  is  very  poor  for  small  n.  For 
example,  if  m  =  3  and  n  =  10,  it  gives  R  whereas  the  exact  value  is  R  =  0.  67. 

It  is  much  closer  to  the  lower  bound  R  >  0.  47  in  Equation  (  22  ).  Even  the  approxima¬ 
tion,  (  23  ),  gives  the  poor  estimate  R  «:i0.4.  Therefore,  it  is  advisable  to  determine 

m 

the  bounds  given  by  Equation  (  21  ) ,  and  to  use  the  exact  expression,  (  20  ) ,  if  they  are  too 
far  apart. 

COROLLARY 

The  term  R^  ^  (0)  -  1,  combined  with  Equation  (  21  ) ,  gives  the  identity 


,  „  n-1  ^  n-1  n-2 

n  =  1+3  — r  +  5  — -  .  — T  + 
n+1  n+1  n+2 


.  (2n-l)  ^  . 


2n-l 


(  26) 


for  all  integer  n. 


This  corresponds  to  the  upper  bound  in  Equation  (  21  ). 
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SECTION  VI 

THE  REDUCTION  FACTOR  AT  THE  PREDICTED  POINT 


The  predicted  point  is  the  first  point  outside  the  given  interval,  x  =  -1  (x  =n  gives 
the  same  results).  The  value  (-1  ,n)  is  derived  in  the  appendix.  Substitution  of 
Equations  (  52  )  and  (  4  )  into  (  9  )  gives 


=  n 

m  n 


1+3 


n-1  (n+1)  ^  5  n-2 

n+1  *  -^2  n+1  '  n+2 

(n-l) 


(n+1)  (n+2) 


(n-l)  (n-2) 


n-l 

n-m 

^n+1 

n+m 

2 

n+1  •  ’ 

n+m 

n-l  ’ ' 

n-m 

(27) 


and 


R  (-1)  =- 
m  n 


,  ^  n+1  _  n+1  n+2 

1+3— T  +  5;^  .  — 5 

n-l  n-l  n-2 


+  (2mtl)  iil 
n-l 


n+m 

n-m 


(28) 


Equation  (  28  )  has  a  remarkable  symmetry  to  (  21  ).  In  a  similar  way,  we  derive  the 
bounds 


(m+l)^ 

n 


<  R  (-1) 

m  ' 


< 


(m+l)^ 

m 


m+1  n+m 
n-l  '  ’ '  n-m 


(29) 


and  the  approximation 


R 

m 


(-1) 


(m+l)^ 

n 


^  m(m+2) 
2n 


(30) 


I 
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SECTION  vn 

THE  ASYMPTOTIC  BEHAVIOR  OF  THE  REDUCTION  FACTOR 

Since  <}>^(x,n)  is  rapidly  increasing  with  x  for  x  >  n,  the  variance  of  the 
extrapolated  polynomial  y(x)  is  rapidly  increasing.  Exact  values  can  be  obtained  by 
lengthy  computations  only.  An  asymptotic  estimate,  however,  is  quite  easy.  For  large  x, 
the  highest  term  becomes  overwhelming  in  <}>  (x,n) ,  therefore 


A  (X)  .  (-1/^) 


(31) 


Similarly,  in  Equation  (  9  )  the  term  with  r  =  m  becomes  overwhelming.  Therefore, 

1  2 


R  (X)  -  X 
m 


2r 


1-  2  . . .  (2m) 


(1-2  ...  m) 


2  '  (n-1)  . . .  (n-m) 


2m+l  n-1 


n 


n-m 


n+1  ’  ’ '  n+m 


(  32) 


or 


R  (x) 
m 


2r 


1>  3  .  ..  (2m-l) 

2-  4  .  . .  2m 


2m+l 


n^-l 


2  ,  2, 

n  -(m  ) 


(  33) 


This  gives  the  asymptotic  estimate 


/.  \  2r  2  2 

TD  /^\  f/  s  n  n 

R^(d)~f(m)  -  ^■••-2-^  ■ 

\  I  n-1  n  -m 


(34) 


The  factor,  f(m) ,  in  Equation  (  33  ) ,  is  bounded  by  2/n  and  3/4.  For  small  m,  it 
is  closer  to  3/4;  for  large  m,  closer  to  2/n.  Often  the  approximation  f(m)  =  0.  7 
will  be  sufficient. 
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APPENDIX 

SOME  PROPERTIES  OF  GRAM-POLYNOMIALS 


A  RECURSION  FORMULA 

For  any  orthogonal  polynomial ,  a  recursion  formula , 


(j.^  (X)  =  a^x(t)^_^  (X)  +  (X)  +  (x) 


(  35) 


exists.  The  coefficients  a^,  and  c^  can  be  derived  by  the  following  arguments: 

we  note  that  (0)  =1  for  all  r;  therefore, 


1  =  b  +  c 
r  r 


(36) 


Multiplication  of  Equation  (  35  )  with  (j>^  ^  (x)  and  summation  over  all  x  gives 


0  =  a 


2 

)  X(|>  (x)  +  b  q 

4=0  ^  r  r-1 


(  37) 


The  same  operation  with  (|)^  ^  gives 


^n-1 
x=0 


(  38) 


Since  <j)  (x)  is  a  symmetric  or  antisymmetric  function  with  respect  to  the  midpoint 
^  2 

(n-1)  /2  of  the  interval,  ^  ^  (x)  is  symmetric  and  it  holds  that 


n-1 

2 


«!•  r  ° 


(39) 
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This  gives 


-1 

=0 


(40) 


for  the  sum  in  Equation  (  37  ). 

To  evaluate  the  sum  in  Equation  (  38  ) ,  we  develop  the  polynomial  x(j>^  ^  M  ^ 
series  of  polynomials; 

x(j»  (x)  =  a(j)  (x)  +  polynomials  of  lower  degree.  (  41  ) 

r  r  “1. 

From  this  it  follows  that 


x*!*  (X)  <1'^_2  (x)  =  aq^_^  .  (  42  ) 

The  term  a  is  obviously  the  quotient  of  the  highest  coefficients  of  x  in  <j)  and 

r 

A  .  ,  which,  from  Equation  (  3  ) ,  can  be  found  to  be 
r-1 

t 

From  Equations  (  36  ) ,  (  37  ) ,  (  38  ) ,  (  40  ) ,  (  42  ) ,  and  (  43  ) ,  we  obtain  a^ ,  b^ ,  and 
c^,  and  have  the  recursion  formula  for  Gram-polynomials: 


(j)^(x,n)  =  - 


2(2r-l) 

r(n-r) 


X 

I 


n-1 

2 


r-1  n+r-1 
r  n-r 


(44) 


THE  VALUE  AT  THE  MIDPOINT 

For  the  midpoint,  x  =  (n-1)  /  2,  Equation  (  44  )  simplifies  to 


^r 


n 

I 


r-1  n+r-1  j  /n-1 
r  n-r  ^r-2  2 


(45) 
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This  allows  recursive  calculation  of  the  even-order  polynomials  for  the  midpoint;  namely, 


2s 


n-1 

2 


2s -1  n+1  n+3  n+2s-l 

2s  *  n-2  *  n-4  *  *  *  n-2s 


(46) 


The  odd-order  polynomials  equal  zero  at  the  midpoint, 

THE  VALUE  AT  THE  PREDICTED  POINT 

We  call  X  =  -1  the  ’’predicted  point.  ”  Since  (-1)^^^  =  (-l)^i !  ,  Equation  (  3  ) 
simplifies  to 


/  =  1  +  +  r(r-l)  ■  (r+1)  (r+2) 

r'  '  l(n-l)  1-  2  •  (n-1)  (n-2) 


r(r-l)  ,  , ,  1  ■  (r+1)  , , ,  (2r) 
1-2  ...  r  •  (n-1) . . .  (n-r) 


This  is  a  hyper  geometric  series,  namely, 


(j)^(-l,n)  =  F(r+1; -r;-n+l;l) 


For  Gauss’  relation, 


[2] 


3 


(48) 


(c-a-b)  F  (a,b;cpc)  -  (c-a)  F  (a-l,b;cpc)  +  b  (1-x)  F(a,b+l;c;x)  =0  ,  (49) 


we  obtain,  for  x  =  0,  since  F  is  finite, 


c-a 


F(a,b;c;l)  =  F(a-l,b;c;l) 


(50) 


And,  if  a  >  0,  by  induction. 


„/  ,  _  c-a  c-a+1  c-1 

c-a-b+1  **'  c-b-1  ’ 


(51) 


3 


See  Formula  6. 


since  F(0,b;c;l)  =1.  Substitution  of  c  =  -n+1,  a  =  r+1  and  b  =  -; 
re-arrangement  gives 


=F(r+l,  -r;-n+l;l)  =  — 


n+2  n+r 

n-1  ' '  ‘  n-r 
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and 
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